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Abstract. In a paper by Maartens, Lesame and Ellis (1998) it was shown that 
irrotational dust solutions with vanishing electric part of the Weyl tensor are subject 
to severe integrability conditions and it was conjectured that the only such solutions 
are FLRW spacetimes. In their analysis the possibility of a cosmological constant A 
was omitted. The conjecture is proved, irrespective as to whether A is zero or not, and 
qualitative differences with the case of vanishing magnetic Weyl curvature are pointed 
out. 
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1. Introduction 

Irrotational dust (ID) spacetimes are perfect fluid solutions of the general relativistic 
field equations characterized by vanishing pressure (p = 0) and vorticity vector {uj a = 0). 
They serve as potential models for the late universe [HE] and gravitational collapse [Sj. 
In an initial- value problem formulation for such models, the Ricci-identity for the fluid 
4-velocity and the Bianchi equations (incorporating the field equations via the Ricci 
tensor) may be covariantly split into propagation and constraint equations [I]. In a 
streamlined setting 5] Maartens showed that these constraint equations are consistent 
with each other (i.e., not overdetermined on an initial hypersurface) and are preserved 
under evolution, generically. However, when further external conditions are imposed 
their consistency is not a priori guaranteed, and rather unexpected conclusions may 
come out. 

In this respect, one of the most promising and natural further restrictions one can 
make is to set the gravito-magnetic tensor H a b equal to zero (a cosmological motivation 
for doing so was given in with the remaining gravito-electric tensor E a b ^ 0. As 
the latter is the relativistic generalization of the tidal tensor in Newtonian theory, 
the corresponding 'purely electric' ID models were called Newtonian-like in One 
might therefore expect to find a rich class of solutions, at least incorporating translates 
of classical Newtonian models. However, in two independent papers |H1 Ej a non- 
terminating chain Ce of consistency conditions for the dynamical variables was deduced, 
which made the authors conjecture that the class was unlikely to extend beyond the 
known spatially homogeneous members and the Szekeres JU] subfamily, which exhausts 
Petrov type D [llj. 
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The natural counterpart to the above is to set E ab equal to zero, with H ab ^ 0. 
These purely magnetic ID models were studied in [7j and were logically called 'anti- 
Newtonian'. Some parallels and differences between Newtonian-like and anti- Newtonian 
universes were pointed out. In both cases no spatial derivatives occur in the remaining 
propagation equations; hence the flowlines emerging from an initial hypersurface evolve 
separately from each other, and solutions may therefore be called 'silent' [JJ, although in 
the literature this terminology is normally used and was introduced ^2] for the purely 
electric case only. Also, the chain Ce for the Newtonian-like case has a direct analogue 
Ch,i in the anti-Newtonian case. However, whereas Ce is identically satisfied in the 
linearized theory (with Friedman-Lemaitre- Robertson- Walker (FLRW) background), 
leading to a linearization instability for Newtonian-like ID models (in the sense that there 
are consistent linearized solutions which are not the limit of any consistent solutions in 
the full, nonlinear theory [8j), it was shown by the primary integrability condition in 
Ch,i that there are no linearized anti-Newtonian ID models at all. On top of this, a 
second chain Ch.2 of integrability conditions exists for the magnetic case, which has no 
analogue in the electric case. Both facts have lead the authors of [7j to conjecture the 
non-existence of purely magnetic ID models, but a proof for it within the full, nonlinear 
theory was not found until now. 

It is usually believed that the presence of a cosmological constant A is of minor or no 
importance to the qualitative outcome of consistency analyses for particular conditions, 
and in this line A was omitted in [HI IHj and [JJ. For Newtonian-like silent models 
however, one of the special cases where one can show that the constraints with A < 
are inconsistent (namely the case where one of the eigenvalues of the expansion tensor 
is 0), surprisingly turns out to admit two new and explicitly constructable families of 
spatially mhomogeneous Petrov type I metrics when A > [15] . This result points out 
that a cosmological constant should preferably be incorporated in consistency analyses. 

In the present paper and in this respect, the generalization to arbitrary cosmological 
constant A of the above Maartens et al conjecture for purely magnetic ID models is 
proved, by showing that the combined conditions of vanishing pressure p, vorticity 
vector u a and gravito-electric tensor E ab are only consistent for FLRW spacetimes. 
The structure of the paper is as follows. In section 2.1 the 1+3 covariant system of 
propagation and constraint equations for general irrotational dust is written down as in 
[7j, with a cosmological constant added. In section 2.2 the origin of the two chains of 
integrability conditions for the purely magnetic case is resumed and it is pointed out how 
a third chain is coming into play. Comparison with the purely electric case is made in 
Section 2.3, and the conjecture for Newtonian-like ID models [El El is properly restated 
in function of the sign of A. The actual proof (in two steps) for the inconsistency of 
anti- Newtonian universes is the content of Section 3. A natural continuation of the 
analysis is pointed out in Section 4. 
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2. Magnetic versus electric ID 

2.1. Basic setting for general ID 

Einstein's field equation with cosmological constant A and perfect fluid source term 
T a b = l^u a Ub +ph a b can be covariantly written as (8ttG = 1 = c): 
1 1 

R ab = ^ (P + 3p)u a ii b + - (p - p)h ab - Ag ab . (1) 

It provides an algebraic definition of the Ricci tensor R ab in terms of the matter 
density p, pressure p and normalized timelike four- velocity field u a (u a u a = —1) 
of the perfect fluid; the corresponding pair of complementary idempotent operators 
(— u a u b , h a b := g a b + u a u b ) projects along and orthogonal to u a , respectively 

The introduction of the following 1+3 covariant projections and operations (see 
EH) considerably compactifies expressions. Herein round (square) brackets denote 
symmetrization (anti-symmetrization) and tr(S') is written for S a a ; for 2-tensors U a b 
and one-forms X a the 'dot-free' notation (SU) ab for S a c U cbl resp. (SX) a for S b X b) is 
used throughout the paper. 

• The spatial part of a tensor S ai '" ar bl ,„ bs is 

S0ai...8r LOl LOr U d\ 1 d B QCl...Cr 1 0\ 

^ 6i...6 s — 11 ci • ■ ■ I 1 Cr^bj • ■ ■ itbs J di...d s - \ Z ) 

A tensor is called spatial iff s S ai '" ar bl ... ba = S ai '" ar bl ... bs . The spatial, symmetric and 
tracefree part S{ a b) of a 2-tensor S ab is 

S {ab ) = s S {ab) - hr( s S)h ab . (3) 

• The spatial permutation tensor is e abc = e [abc] = r] abcd u d , where i] abcd = r] [abcd] is the 
spacetime Levi-Civita permutation tensor. One then defines for 2-tensors S ab , U ab : 

[S,U} a = e abc (SU) bc . (4) 

When S a b and U ab are spatial and symmetric, one can easily show the identity (see 
also (A3) o/0) 

[S 2 , U] a = 2 [S, SU] a = tr{S) [S, U] a - (S [S, U)) a , (5) 

with SUab = (SU)(ab). In this case [S, U] a = is equivalent to (SU) ab being 
symmetric, which is still equivalent to (SU — US) ab = 0, whence S a b and U a b have 
a common eigenframe. 

• The electric part E ab and magnetic part H ab of the Weyl tensor C abc d are 

E a b = Ei^ab) = C acb dU c u d 1 H ab = H(a b ) = e c d a C cd be u e . (6) 

• Let V be the spacetime covariant derivative attached to the unique metric 
connection. The (covariant) time derivative T and the (covariant) spatial derivative 
D are defined by: 

T(S) a ^ bi ba =u c V c S"r bl ... ba , (7) 
P, S'" ,„..,,_ = s V c ^---<V.. bs . (8) 
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Normally the notation S ai --- ar bl ...b s is used instead of T(S) ai - ar b h , but the latter 
emphasizes on the action of T as a (degree derivative) operator and this is more 
convenient for what follows. Definition (JHJ) tells that D maps each (r, s)-tensor 
S ai '" ar b 1 ...b s to the spatial part of the (r, s + l)-tensor V c S ai '" ar b 1 ...b s - It gives rise to 
the covariant spatial divergence (div) and curl operator, acting on vectors V a and 
2-tensors S a b as: 

div V := D a V a , curl V a := e abc D b V c (9) 
div S a := D b S a b, curl S aft := e cd(a D c S b) d ; (10) 
For spatial 2-tensors S a b, U a b one has the identity 

div[S,U] = tr([/curlS) -tv(ScmlU). (11) 
• The splitting of Vm into different kinematical quantities [3] : 

U a -b = D b U a - U a U b 

■■= o-ab + \^h ab + e abc u c - u a u b . (12) 

Herein a a b = u^ a -b) is the shear tensor, 6 = u a - a the volume-expansion scalar, 
uj a = \e ahc Ub-c the vorticity vector, and u a the acceleration vector, which are all 
spatial. 

For this introduction, the abstract index notation {H] was most transparent, but 
an index-free notation turns out to be practical for calculation purposes. In particular, 
E a b, H a b,o- a b,h a b,u a ,u a are denoted E, H,a,u,u from now on. Moreover, since all 
appearing 1-forms and 2-tensors will be spatial the symbol h will be dropped everywhere. 
The notation for S/ab\ becomes S (cf. [18 ), but when S is symmetric this is written as 
S — ^trS for computational convenience. In this notation, the M-linear operators T and 
D satisfy in particular the Leibniz rules 

T{XY) = T(X) Y + X T(Y), D{fg) = (Df) g + f Dg (13) 

for 2-tensors or scalar functions X, 2-tensors, 1-forms or scalar functions Y, and scalar 
functions f,g. In the present paper we deal with non-vacuum irrotational dust (ID), 
i.e. solutions of for which 

p = 0, p^O, u = 0, ii = 0, (14) 

where u = is a mere consequence of p = 0, p ^ and the second contracted Bianchi 
identity. With u = one has 

T[S,U] = [T(S),U] + [S,T(U)]. (15) 

A commutator rule for T and D, written in appropriate form and for the action on 
scalar functions /, is also needed. With u — it = it reads [3]: 

T(Df) = D(T(f))-(a + h)Df. (16) 

For perfect fluids in general the Ricci identity for u and the Bianchi-identities, 
hereby substituting (JTJ) for the Ricci tensor, are basic equations to be satisfied. They 
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may be covariantly split into a set of propagation equations (involving the covariant 
time derivative) and a set of constraint equations jUE|. In the case of non- vacuum ID 
(JUJ) the sets of non-identical equations read: 

• Propagation equations: 

T(p) = -pO, (17) 
T(#) = -^ 2 -tr(a 2 )-ip + A, (18) 

T(a) = -ha - a 2 + hr(a 2 ) - E (19) 

T(E) - curl// = -BE + Za~E - ^pa, (20) 
T{H) + curlE = -6H + 3aH (21) 

• Constraint equations: 

2 

C 1 = diva - -DO = 0, (22) 
3 

C 2 = curia - H = 0, (23) 

C 3 = div^ - -Dp - [a, H] = 0, (24) 
3 

C 4 = divH + [a, E] = 0. (25) 



On considering (JTHJ) and (J2BI), fl2H, (122) and (j23), one sees that within the class of 
non-vacuum ID spacetimes, the FLRW case 

Dp = D9 = 0, a = E = H = 0. (26) 

is covariantly characterized by vanishing shear tensor (cr = 0) or, equivalently, by Petrov 
type (E = H = 0). 

In jH] it was shown that the time derivatives of each of the C l ,i = 1..4, are 
differential combinations of the C J themselves, such that (}2*2*j) -()25 |) evolves consistently 
and doesn't give rise to new constraints, generically. It was also shown that actually, C 4 
is a differential combination of C 1 and C 2 : 

C 4 = ^curlC 1 -divC 2 . 

However, if extra covariant conditions are imposed, new constraints may be generated 
and their consistency should be investigated. 



2.2. Purely magnetic ID 

In this respect, the constraint E = was studied as the natural counterpart for H = 
in [Jj. The remaining variables are the scalars p, 9 and the tensorial quantities a,H. 
As put forward in the introduction, one sees that no spatial derivatives occur in the 
propagation equations (fl"Tj) - (fT9"j) and (|2T|). such that the purely magnetic ID models 
are 'silent'. Actually, with E = 0, (fT7|) - (fr9^) forms a closed subsystem determining the 
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evolution of the matter variables p, 9 and a, from which the evolution equation (|21j) 
for the gravito-magnetic field H is decoupled. The latter equation may also be seen as 
arising from the curl of the algebraic constraint E = put on the whole of the system 
(fT7j) - (j23j) . In the same viewpoint and by virtue of (}2T)j) and (J2H), the time derivative 
and the divergence of E — gives rise to two independent constraints which are not 
preserved under evolution: 

C 5 = curl/7 - l -po = 0, (27) 

C W = -Idp-[ct,H} = 0. (28) 

It was not noted in [2j that, although the sum equation C 3 = of div E = and 
q3,h = q j s consistent under evolution, this is not necessarily the case for both equations 
separately. As E = has no influence on C 1 and C 2 , and as (f2Tj) is an identity, the 
constraints (|2^j). (J2llj) and (J23j) do evolve consistently as in the generic case. 

The time derivative of (|27jl . modulo (|22|) and (|27j) . yields a primary integrability 
condition £ = whose repeated time differentiation leads to a non-terminating chain of 
new constraints Ch,i after . The linearization of £ around an FLRW background (|2*Bj) 
reads \p9a, such that linearized non-vacuum ID models with E = are either FLRW 
or satisfy 9 = 0,p = 2A, as is seen from ()18|). Hence in the case A = no linearized 
anti-Newtonian universes exist at all [7j. Also note that, in the light of the question 
whether to incorporate a cosmological constant or not, the same conclusion cannot be 
drawn for general A). However, for the purpose of proving inconsistency in the full 
non-linear theory, the chain Ch,i is even more unmanageable than its purely electric 
analogue Ce (see Section 2.3). 

Analogously, for the time derivative of C 3,H one finds after a short calculation, on 
using (HHJ), (USD, (HH, (j2U) and ©: 

T(C^) = Qa - jj*) C 3 <" + \j, 

J = p D9 + (\a -V) Dp. (29) 



2 3 , 

Henceforth, J = 0. As E = 0, and C 3,H = is the translation of div E = whereas 
C 5 = is coming from E = 0, the same condition should be found on calculating divC 5 
(see also (A13) in 0). Indeed, in [7j it was calculated that J 

div C 5 = -curl C 4 - (a + -9) C 3 - -C 1 - -J. (30) 
2 3 2 3 

The successive time derivatives of := J = lead to a second chain of constraints 

C H)2 = (J"W = 0)^ , with := T l (:7). In [7j and J^ 2 ) were explicitly calculated 

(for A = 0) and the general structure of J^ 1 was given, but apparently the consequence 

of this chain was overlooked. In Proposition 1 of Section 3, it will be shown (for arbitrary 

A) that Ch,2 implies Dp = 0. But this gives rise to a third chain C# 3 of constraints, as 

J In 0, the — aC 3 term was forgotten in the right hand side. 
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follows. With Dp = 0, (|28|) reads [a, H] = 0. Taking the divergence of this and using 
(fTT |) . (Pjl and (|27jl . one finds: 



and hence C H ,3 = (£ W = 0)£ , with £W := T i (/C)§. In Theorem 1 of Section 3 it is 
finally shown that within the class of ID's with E — 0, Ch,3 can only be satisfied for 
a = H = 0, i.e., by the FLRW subclass. 

2.3. Comparison with purely electric ID 

Analogous to the purely magnetic case, (jl7|) - (pnj) with H = forms a 'silent' system 
of evolution equations for the remaining variables p, 9, a, E, but here (J2U|) is coupled to 
P7|) - (Pll via the E-term in (JTHJ) - Now H = (together with F = 0) and div H = 
translate into 



C 1 is independent of H, and C 2 , C 3 contain H only algebraically; hence the corresponding 
constraints evolve consistently. In analogy with the E = case, the time derivative of 
(HUD), modulo the constraints (jHJ)-(j2U) and (jH2|), gives rise to a new tensorial condition 
7i = and an indefinite chain Ce = (T l (H,) = 0)^ of constraints, identically satisfied 
for Petrov type D and spatially homogeneous models, but not in general. However, 
in contrast to the magnetic case, the linearizations of the T l {H) around an FLRW 
background are identically zero; hence Newtonian-like silent models are subject to a 
linearization instability, cf. Introduction. The difference is essentially due to the extra 
term ^pcr in ()27|) within the magnetic case, not present in ()32|). The same absence 
makes that the divergence of C 5 doesn't lead to a new constraint in the electric case: 



By the analogy of the remark made for the magnetic case, the propagation of the 
constraint should not lead to a new condition either. Indeed, this can easily be 
checked on using ((T3J), (JTHJ), (j2~Uj) with H = 0, and Moreover, it was checked in [TE] 
that even for general ID the condition div if = 0, or equivalently [a, E] = 0, is consistent 
under evolution, this making use of the time derivative of (fT^j) instead of using (|2*U|) ||. 
Hence there is no analogue here for the chain Ch,2 obtained in the magnetic case, and 
the difference may also be traced back to the fact that, as in electro-magnetism, Dp 

§ Notice that, without the extra fact Dp = 0, the divergence of divC 3,H already incorporates a 
Laplacian D 2 (p) 1 which makes the direct chain T l (divC 3 H ) — rather unmanageable. 
|| Essentially, this comes down to showing [a, T(a)] = [a, T 2 (a)\ = 0. The first equation follows trivially 
from 1)19(1 and [cr, E] = 0; the second one then immediately follows from time differentiation of the first 
and applying (15), but strangely enough the authors of 16 turned to a much longer but equivalent 
reasoning in an orthonormal tetrad approach to show this. 



/C = tr(iJ 2 )--ptr(a 2 ) = 0, 



(31) 



C 5 = curl E = 0, 
C*> E = [a, E] = 0. 



(32) 
(33) 




(34) 
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influences the divergence of the electric, but not the magnetic field. Hence div H = 
does not lead to direct information involving Dp. 

Finally notice that a counterpart for the third chain Ch,3 doesn't exist either: the 
divergence of leads to an identity, since curl E = and since (|2*3*|) now reads 
curia = 0! 

Since [a, E] = 0, (fTTjl - ff^Uj) gives rise to a autonomous dynamical system S in p, 9, 
and two Weyl and shear eigenvalues (or combinations thereof). Analyses for A = in 
an orthonormal Weyl eigentetrad approach [H| or coordinate approach ;9] both led to 
the same conclusion, which remains valid for general A: for the spatially inhomogeneous 
Petrov type I subclass C/ iS i of the purely electric silent models, the chain Ce eventually 
translates into one or more indefinite chains of polynomial constraints on S. In |13j . 
four (compact) relations on the variables x of S were found, by which all these chains 
become identically satisfied. Hereby A seems to play the role of a bifurcation parameter 
for S: for A < and real x, these relations imply Petrov type (FLRW), while for A > 
members of C/, s i emerge. Further analysis by computer (details of which will be given 
elsewhere) strongly indicates that these are the only members of Cj )S i for A > 0, whereas 
there are probably none for A < 0. This is a generalized statement of the conjecture in 
[HI HI] • Still no definitive proof has been found, mainly because of the massiveness of the 
polynomials. 

3. Anti-Newtonian universes do not exist 

As outlined in Section 2.2, it is shown in this section that non-vacuum irrotational dust 
models with E = imply Dp = (Proposition 1) and hence are FLRW (Theorem 1), 
i.e., that anti-Newtonian universes do not exist. Hereby, all statements and expressions 
are implicitly assumed to be related to, resp. live on, a small open subset U of a 
spacetime. Since the whole of the reasoning only involves a finite number of polynomial 
combinations Fi of tensorial quantities on U, one can always assume U small enough 
such that for all couples (F^Fj) either i*$(p) = Fj(p),\/p G U (denoted by Fj = Fj) or 
Fi(p) ^ Fj(p),Vp e U (denoted by F t ^ Fj). 

It should be remarked here that purely magnetic irrotational vacua (E = 0,uj = 
0,p = p = 0, with or without cosmological constant) were shown to be inconsistent, first 
by Van den Bergh jTJj in an orthonormal tetrad approach, and later in a more transpar- 
ent 1+3 covariant deduction (thereby generalizing the result from vacua to spacetimes 
with vanishing Cotton tensor) by Ferrando and Saez ^H]- Hence the vacuum case p = 
may be excluded from the subsequent analysis. 

Some additional preparations should still be made. For the first step of the proof 
of Proposition 1, = T'\J) have to be calculated up to % = 6. Denote x fe for tra^ 
(k G N). By (|Tnj) with E — 0, T(er) is an element of the commutative (!) polynomial 
algebra (over spacetime functions) generated by a itself, such that T(a k ) = ko- k ~ 1 T(a). 
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Taking the trace and substituting for T(a) one finds 

T(x k ) = k Xk ~ x k+i + ^x 2 x k -^j . (35) 

Now the key remark is that all of the x k for k > 4 are polynomially dependent of x 2 
and x 3 (according to the formulation in [7j w.r.t. these traces, this seems to be exactly 
what was overlooked). Indeed, on using Newton's identities up to power three, applied 
on the eigenvalues of a, the Cayley-Hamilton theorem for a reads 
.11 

cr = -x 2 a + -x 3 . (36) 

On multiplying (jHBj) with a k and taking the trace (or just taking Newton's identities for 
powers higher than three) one gets 

x k+3 = ~x 2 x k+1 + ~x 3 x k , fc=l,2, .... (37) 

Since X\ = tr a = this reads 24 = \x 2 2 for k = 1. Hence for = 2, 3 (|33|) becomes: 

T(x 2 ) = - -6x 2 - 2x 3 (38) 

T(x 3 ) = - 26x 3 - X -x\. (39) 

Hence (|3*8 )) -(|3'9" j) together with (fTTj) - (P"9"j) forms an autonomous dynamical system Sh 
in the scalar invariants p,9,x 2 ,x 3 . Starting from ()29jl . it is easily seen that for the 
computation of any one only needs the basic operations (fTSj) and (fTB^l. together 
with (fTT^) and the equations of Sh as substitution rules for T(a),T(p),T(9),T(x 2 ) and 
T(x3). All w ill be elements of the module generated by Dp, pD9, pDx 2 and pDx 3 
over the (commutative) polynomial ring R = Q[p, 9, x 2 , x 3 , A][s]. In particular, one finds 
for JO) and J<®\ 

= - pDx2 _ (l a + 5^ pZ ^ _ Q (p + A) + a 2 + 2 _ da _ I _ Dpj 
J® = 2 pD{x 3 ) + + P D{x 2 ) + (lp - ^A + 4a 2 + ^9a + 2x 2 + ™9 2 ^J P D{9) 

+ ((^ + 1 ^)P + * ~ r) A + 3 ^ 3 + 4 ^ 2 " (y " - - ^ - y^ 3 ) ^(P)- 

Hence J = = = can be solved for pD9, pDx 2 and pDx 3 , which expresses 
them as elements of the module generated by Dp over R % 

pD9 = (le - \a\ Dp, (40) 
pDx 2 = - (\ 2 + l -6o - X -x 2 + ~(p + A)] Dp (41) 



(^ + ^-^ 2 + i(p + A))D P , 



(42) 



142(1 and 144(1 were added for comparison with the corresponding formulas in the case A = from 
0. Apparently there was a typo regarding the coefficient of x-ioDp in 1(44(1 . 
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pDx 3 = ( A(ia + l -B) - ^ - Ua 2 + l -x 2 a - ^x 2 9 + \x 3 ) Dp (43) 

A(-a + -0) - —a 3 - -6a 2 + -x 2 a - -x 2 9 + —x 3 ) Dp. (44) 
v 4 6 ; 16 8 8 8 16 1 ' 

The result of Proposition 1 is that anti-Newtonian universes necessarily satisfy 
Dp = 0, and hence [a, H] = from (|28|). As a and H are both spatial and symmetric, 
this implies that (a) o~H is symmetric, such that one can write (|21|) as 

T(H) = -6H + 3aH -tr(aH), (45) 

and (b) that a and H commute as (l,l)-tensors. Hence T(o~) and T(H) are elements of 
the commutative polynomial algebra (over spacetime functions) generated by a and H, 
such that T(a l Hi) = a l ~ 1 H^~ 1 (iHT(a) +j aT(H)) for arbitrary i and j. Substituting 
for T(a),T(H) by (0,(113) and taking the trace yields: 

2i i 
T ( x i,j) = ~ ( y + jW x i,j + (3j - + g^2,o^i-i,j - (46) 

wherein Xjj is written for tr(<7*-£P). On multiplying (}3l)j) with a'if J and taking the trace 
these are found to be related by 

1 1 

x i+s,j = 2 X 2,o^+ij + ^XsjaXij, i,j = 0, 1, 2, . . . . (47) 

For the first step of the proof of Theorem 1, K^' = T l (/C) needs to be computed up 
to i — 5. Starting from /C = x ,2 — ^px 2 ,o, inspection of (j4T)j) for j = 2 learns that 
KM\i = 0..5 will be of the form r s ® i y %i,2+ other terms. Hence KM' = determines 
Xi 2 as a polynomial in p, 9, A, Xj t2 (j < i) and some of the Xk,o,Xi t i. Substitution in a 
particular order of identities (|4Tj) and obtained expressions yields compact equations 
u[i] — 0, i — .. 5, where 

1 . , 

u = x 0:2 - -px 2fi (48) 

1 

ui = -x 2fl p6 + 6xi j2 + px 3 , (49) 
6 

U 2 = - ~X3,OP0 + 7;PX2,0 2 - T^X 2fi p 2 + ~X 2 . pA + 30X2,2 - 12X1,! 2 (50) 

3 3 12 6 

u 3 = -p 2 x 3 , + -x 2fi Ap9 - 2x 2}0 2 p9 + 6x 3 , px 2 ,o - 108x1,1X2,1 - x 3 , pA(51) 

3 4 
u A = - 18xi,i 2 x 2 , - x 2fi 3 p + -x 2 , V - 12px 3 , 2 - -x 3fi Ap9 + 6x 2fi p9x 3fi 

2 115 

+ -Ap# 2 x 2 , - ~^x 2fi p 2 A + -x 2 , A 2 p - -x 2 , 2 pA - 216x 2 ,i 2 (52) 

29 1 2 

w 5 = 36xi,i 2 x 3 , - —x 2 /Ape - -x 2fi Ap 2 9 - -x 2 /p6 - 36x 2 ,o£i,ix 2 ,i 
10 6 6 

- 8x 30 2 p9 + 6px 2 , 2 a; 3 ,o - ^-a^o^op 2 + 7P 2 Ax 3 , - |x 3 , A 2 p 

2 6 3 

+ ^Ax 2 , # 3 p - -^A0 2 px 3 , + l5x 2fi x 3fiP A + -x 2>0 A 2 pd (53) 
27 9 6 
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Remark. The calculation of J~^\i = 1 .. 6 and u[i], i — 0.. 5 were preformed in the 
standard commutative environment of Maple (version 9.5), as allowed by the index-free 
notation and in the present situation. The used Maple code which reflects how the u[i] 
were obtained, is explicitly given in Appendix A. 

The second step of the proofs of both Proposition 1 and Theorem 1 require two 
independent combinations of the non-zero shear eigenvalues a k and of the magnetic 
Weyl eigenvalues H k each (k = 1, 2, 3). The following choice was made: 

3 1 

si :=--(?!, h := -(cr 2 -<7 3 ), (54) 

h-^—Hi h:=^{H 2 -H 3 ), (55) 
with inverse relations 

2 1 1 

Cl = --Sl, CT 2 = -S 1 +t 1 , 03 = -S]. - ti, (56) 

H 1 = ~h 1 , H 2 = ^h 1 + k 1 , H 3 = ^h 1 -k 1 . (57) 

Lemma 1. For a purely magnetic ID model (a) si = implies FLRW and (b) 
t\ = csi, with c a constant function, implies that the shear tensor is degenerate. 

Proof. The shear tensor is a diagonalizable (1,1) tensor. Let P& be the projector 
on any eigenvectorfield corresponding to the eigenvalue a k . Then it is straightforward 
to deduce that 

T(a k ) = tr(P fe T(a)) = ~9<r k -a 2 k + hi(a 2 ). (58) 

where the second equation is valid for the specific situation (|19|) . From this it follows 
that 

T{ Sl )= - \es x + ^-tl (59) 

T{h) = - ^(sx + 9). (60) 

(a) If S\ = then also t\ = by (|59|). such that the ID model is shear- free and 
hence FLRW. (b) The time derivative of t\ — cs\ = 0, substituting t\ for csi, yields 
c(c— l)(c+ l)sf = 0. The case si = yields FLRW, whereas c = or c = ±1 is exactly 
saying that the shear tensor is degenerate, according to □ 

Finally, the following notations and properties concerning multivariate polynomials 
are used. Consider the polynomial ring Q[x], where x stands for the n-tuple (xi, . . . , x n ). 
For two elements pi and p 2 of Q[x], gcd(pi,p 2 ) is written for one of the two normal- 
ized forms of their greatest common divisor (i.e., for which all coefficients are relatively 
prime integers); hereby gcd(p\,p 2 ) ^ 1, resp. gcd(pi,p2) = 1, indicates that p\ and 
p 2 have, resp. do not have, a common factor. Write Kes(pi,p 2 ; Xi) for the resultant 
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of pi and P2 w.r.t. the variable Xj. It is well known that Res(pi,p2',Xi) is a polyno- 
mial over Q in the variables Xi, . . . , Xj_i, Xj + i, . . . , x n , which is an element of the ideal 
of Q[x\ generated by p\ and P2, i.e., there exist polynomials k±,k2 E Q[x] for which 
Res(pi,p 2 ; Xj) = k x p\ + k 2 p2- Hence, if x° is a root of p\ and p 2 (i.e., a solution of 
Pi(x) = P2(x) = 0) then (x°, . . . , x°_ l5 x° +1 , . . . , x°) is a root of Res(pi,p 2 ; Xi) (see e.g. 
[T§]). The following properties are perhaps less known but sometimes help to avoid 
cumbersome computations at the end of a proof by resultants: 

Lemma 2. (a) For pi,p 2 E Q[x] and any Xj, Res(px,P2j %i) is the zero polynomial 
if and only if gcd(pi,p2) ^ 1- 

(b) When n = 2, and if Pi(xi, x 2 ) and p 2 (xi, x 2 ) are homogeneous polynomials of their 
arguments, then the system {pi = 0,p 2 = 0} 

• only has the trivial solution (xi,x 2 ) = (0,0) when gcd(pi,p 2 ) = 1; 

• is equivalent with the equation gcd(px,p 2 ) = when gcd(pi,p 2 ) 7^ 1- 

Proof, (a) see [19J, p 158; (b) is an almost direct consequence of (a). □ 

With these preparations the main result may now be proved in two steps, as indi- 
cated. 

Proposition 1. An ID spacetime with E = satisfies Dp = 0. 

Proof. Calculating . . . , substituting for D6, Dx 2 and Dx 3 by flUJ), (JUJ) 

and one gets four respective equations pi Dp = with pi E R = Q[p, 9, x 2 , X3, A] [a], 

i = 1 .. 4, i.e., Dp is in the kernel of the (l,l)-tensor p^. In writing out these equations 

in their components w.r.t. a shear eigenframe, one gets twelve equations pijDjp = 0, 

i — 1 .. 4, j — 1 .. 3 (no summation over j), where pij = Pij{p, 0, Si,ti, A) originates from 

th ... 
Pi by substituting the j shear eigenvalue o~j for cr, o\ + cr 2 l + o~z for Xj and finally the 

right hand sides of (foljj) for c^, <r 2 , cr 3 . 

Suppose Dp 7^ 0. Then there exists a non-vanishing component Djp. Whitout loss 

of generality we may assume Dip 7^ 0. Then pn = for i — 1 .. 4. One has e.g.: 

P11 = -t^A 2 + (^ Sl 2 + 2t, 2 + + ^ - V)A + (ls x 2 + ^ 2 )p 

-4si 2 ti 2 - 4si#ti 2 - 3ti 4 
P21 = ((~B + ^i)p + ^ - 20^ 2 - y^! 2 - fe 2 Sl + ^9 3 - 4^ Sl )A 

,20, 4 , a2 . 35 „ 2 10 o 245, 2 26 2 . 

+(— 9 - -si A 2 + Sl 2 + — si 3 et x 2 + — ti 2 si p 

v 9 9 ' V 9 9 18 9 JH 

+ie tl % - - Sl % 2 + — e 2 s l t l 2 + — Sl V + 40V. 

3 3 3 

Now compute the resultants qi(9, Si, ti 2 , A) := Res(pn,pji; p), z = 2. .4 and next the 
resultants r^si, t 2 , A) := Res(g 2 , Qi] 0), i = 3, 4. One gets equations of the form 

r 3 ( Sl , t, 2 , A) ee A 2 (3t 1 2 - A) 4 (2 Sl 2 + It, 2 + A) 4 r 3 *( Sl 2 , t, 2 , A) = 0, (61) 
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r 4 ( Sl , t x \ A) = A 2 (3tx 2 - A) 4 (2 Sl 2 + It, 2 + A) 4 Sl r:( Sl 2 , t x \ A) = 0, (62) 

with gcd(r^,r^) = 1. If vanished then the spacetime would be FLRW according to 
Lemma 1 (a), hence Dp = 0, which is a contradiction. Thus we still have to investigate 
the following four cases: 

(i) A = . One has si, £i, 0) = t 2 cii(#, si, ti) = 0,2 = 2. .4, with gcd(aj,Oj) = 1 
for i ^ j. For Pj(si,ti 2 ) = Res(a 2 , a^; #), i = 3,4 one computes that gcd(P 3 ,P 4 ) = 
^i 2 (7ti 2 +2s 1 2 ) 2 (35ti 4 +42t 1 2 s 1 2 +16si 4 ). If t x ^ then s x and tx must be solutions of 
gcd(P 3 , P 4 ) = because of Lemma 2 (b); if ti = thenp n (p, 9, Si, 0, 0) = \psx 2 = 0. 
Thus both possibilities imply a = 0, which leads to the contradiction Dp = 0, as 
above. 

(ii) A = 3ti 2 . One has fc(0, Si,ti, 3t x 2 ) = + Si) 2 &i(#, si, *i) = 0,i = 2. .4, with 
gcd(6j, &j) = 1 for z 7^ j. If ti = then A = 0, impossible since (i) has just 
been excluded. One has pu(p, — si, si, ti, 3£x 2 ) = fp(si 2 + 5ti 2 ) = 0, such that 
# + si = again leads to a = and contradiction. If £i 7^ 7^ 6> + si then 
bi(9,si,ti) — 0, i = 2. .4. But for Qi(si,£i) = Res(&2, #), z = 3,4 one computes 
that gcd(Q 3 ,Q/i) = 1. Hence a = by Lemma 2 (b), contradiction. 

(iii) A = — 2si 2 — 7ti 2 . This case is exactly equivalent to the vanishing of the leading 
coefficient of pxx, seen as a polynomial of p, i.e., one has Pu(p, 9, Sx, tx, — 2s% 2 — 
7ti 2 ) := F(9, si, ti 2 ), the exact form of which is compact but not relevant. Set 
Ri(p, 9, si, ti 2 ) = pa(p, 9, s\, tx, — 2si 2 — 7ti 2 ), i = 2. .4, compute the resultants 
Ui(9, Sx,tx) = Res(P 2 ,P;;p) and then the resultants Vi(sx,tx) = Res(P, Iff, 9), 
i = 3,4. One finds gcd(V 3 ,V 4 ) = (It 2 + 2s 2 ) 3 (s 2 + 5t 2 ) 4 , such that a = via 
Lemma 2 (b), contradiction. 

(iv) r% = r 4 = 0. Then, according to Lemma 2 (a), si 2 and t\ 2 are necessarily solutions 
of the non-trivial homogeneous polynomial equation Res(rg,r 4 ; A) = 0. Hence the 
ratio {tx/ Sx) 2 must be constant (the case s\ = is excluded), which can only be 1 
or by Lemma 1 (b). Now one computes that gcd(rg(si 2 , Sx 2 , A), r 4 (si 2 , Sx 2 , A)) = 
s 2 (3si 2 - A) and gcd(rg(si 2 , 0, A), r^si 2 , 0, A)) = s 2 A 3 , which should be zero 
according to Lemma 2 (b); but this is impossible as (i) and (ii) have already been 
excluded. 

We conclude that the assumption Dp 7^ was false and this finishes the proof. □ 
Theorem 1. Any ID spacetime with E = is FLRW. 

Proof. From Proposition 1 and C A ' H = one has [a, H] = and hence the 
chain C H ,3 = (£ w ee T*(/C) = 0)£ , with K given by (jSH). From the above, the 
Ui,i = 0..5 given by (P5j) -(|53 |) must vanish. As x 2 ,o = is equivalent with a = 0, 
one may divide out factors which are powers of 0:2,0 (or p) in equations F = 0. 
For any polynomial p write p* for the polynomial p from which all such powers are 
taken out, and write x for the variable set (^2,2, £3,05 x 2j x, ^1,2, ^2,0? ^0,2)- One 
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subsequently computes 0,(x, p, #) = Res*(w 2 , uf, A), Pj(x, p) = Res* (ui, Of, 9) and 
Qt(x) = Res*(wo, Pi] p) where i runs from 3 to 5. Now in the Qi perform the substitutions 
= &i l Hi 3 + a 2 l H 2 3 + o~3 l H3 3 for the index couples {i, j) appearing in x, followed by 
the substitutions (jo^j) -(foT j) . This yields polynomials Ri(si,ti, hi, An), i = 3 .. 5. The case 
Si = leads to FLRW by Lemma 1 (a), so we may assume s\ 7^ and scale t\ with s\ and 
hi, ki with s 2 to ease the symbolic computations. This boils down to the substitution 
si = 1 in the Now for Y^ti 2 , hi 2 ) := Res(i?i(l, £1, hi, ki), Ri(l, £1, hi, hi); An), % — 2, 3 
one computes that Z : = gcd(Fi,F 2 ) = (3£i 6 +75£i 4 -15£i 2 +l)(3£i 2 +l) 29 /ii 24 . The second 
factor is positive definite, and so is the first, as 75£i 4 — 15£i 2 + 1 = 75((£i 2 — j^) 2 + 3^). 
This leaves two cases: 

(i) hi = . One has Ri(l, t h 0, h) = k^R^ti 2 , h 2 ) a nd R 2 (l, £1, 0, An) = 
kiR* 2 (ti 2 ,k 2 ), where gco\(R{,R* 2 ) = 1. Thus Res(R*, R* 2 ;k 2 ) is a non-trivial 
polynomial in the scaled variable £i 2 by Lemma 2 (a). Hence, if ki 7^ then 
{ti/si) 2 must be a constant, which can only be 1 or by Lemma 1 (b). But 
one has Rl(l,l,k x 2 ) oc ^^(-S^ 2 + 528), 1, k x 2 ) oc ki 6 (-31ki 2 + 2544), 
and Rl(l,0,ki 2 ) oc ki{-k 2 + 3), R* 2 (l,0,ki 2 ) oc ^(k 2 + 3), which leads to 
contradiction. If An = then H = and the spacetime is FLRW. 

(ii) hi 7^ 0. Set Zi := Fj/Z, i = 1, 2. Analogously as for /i! = 0, Res(Z 1; Z 2 ; hi) is a non- 
trivial polynomial in the scaled variable £i 2 by Lemma 2(a), and hence (£i/si) 2 must 
be 1 or by Lemma 1 (b). But one computes that gcd(Zi(l, hi 2 ), Z 2 (l, hi 2 )) = 
gcd(Zi(0,hi 2 ),Z 2 (0,hi 2 )) = 1, i.e. Zi(l,h x 2 ) and Z 2 {\,hi 2 ), resp. Zi(0,hi 2 ) and 
Z 2 (0, hi 2 ), have no common solution, which leads to contradiction. 

All possible cases lead to FLRW or a contradiction, and this finishes the proof. □ 

Remark. The reason why the analysis was first performed in scalar invariants, and 
not in the variables Sx,£i, hi, An from the start, is twofold. Firstly, in Proposition 1 the 
use of D(x 2 ), D(xs) instead of -D(si), D(ti) rules out the special case of degenerate shear 
(the transformation between both sets has a Vandermonde-like determinant). Secondly, 
and w.r.t. the chain Ch,3, the compactness and low degree of the polynomials Ui, i — .. 5 
is striking: when expressed in si,ti,hi, An the number of terms increases by a factor three 
to four. Although eventually this is of no importance here, this may be kept in mind for 
other consistency analyses, when one has to make a choice for variables to be eliminated. 

4. Conclusion and discussion 

It was proved in the full, nonlinear theory of general relativity that purely magnetic 
irrotational dust (ID) models ('anti-Newtonian universes') do not exist, irrespective as 
to whether a cosmological constant A is present or not. This was conjectured in [Z| for 
A = 0. The proof was mainly done in a 1+3 covariant approach, not making use of the 
Gauss equation for the 3-Ricci curvature of 3-surfaces orthogonal to u (which yields extra 
information within a tetrad formalism). The analysis for the purely magnetic models was 
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repeated for the purely electric ones, and this summarized the main differences between 
both cases, as partially pointed out in The conjecture for spatially inhomogeneous 
Petrov type I Newtonian-like ID models from [HI E] was properly restated in function of 
the sign of A. 

As a cosmological constant can be interpreted as a negative constant pressure, ID 
spacetimes are related to the class of irrotational models with an equation of state. In 
this respect, purely electric or magnetic models for which the flowlines form a geodesic 
congruence are most closely related to the analysis here. This will be dealt with in a 
forthcoming paper. 
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Appendix A. 

The following Maple code generates the expressions P%j) - (j53j) : 

x[4,0] := l/2*x[2,0] 2 : 
for i from 3 to 5 do 

x[i, 1] := 1/2 * x[2, 0] * x[i - 2, 1] + 1/3 * x[3, 0] * x[i - 3, 1]: 

CHx2[i] := x[i, 2] = 1/2 * x[2, 0] * x[i - 2, 2] + 1/3 * x[3, 0] * x[i - 3, 2] 

od: 

K[0] := x[0, 2} - \ * p * x[2, 0] : u[0] := K[0}: 
verg:=x[0,2] = |*p*x[2,0]: 
for i to 5 do 

K[i) :=factor(T( J fsT[z - 1])): 
v[i] :=factor(subs(verg,i^[z])): 
vgl[i]:=isolate(u[i], x[i, 2]): 

if i >= 3 then u[i]:=fa.ctor(svLbs(CHx2[i],veYg,K[i})): 

else u[i] := v[i] 

fi: 

verg:=verg union {vgl[i]}: 
od: 

Herein, T is a differential operator which acts on objects x^j according to (j46j). which 
obeys (fT7|l - (fn?]l with E = and trcr 2 = x 2 ,o, and (g^J with tr(aH) = x 1A . 
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